A Cost-effective Ilu Preconditioner for Weather Simulation
نویسندگان
چکیده
To date, the most efficient solver used in the weather sciences for the resolution of linear system in numerical weather prediction is the generalized minimal residual method called GMRES. However, difficulties still appear in matrix resolution when the GMRES iterative method is used without an appropriate preconditioner. For improving the computation speed in numerically solving weather equations, we investigate the use of Incomplete LU(ILU) factorization preconditioner. This is a fast algorithm for solving large scale linear equations. Some strategies for choosing the level of fill and threshold are described in this paper. As an example, we use the GMRES iterative algorithm to solve the finite difference equations of shallow equations and analyze the results that are obtained in preconditioned and un-preconditioned, respectively. It is shown that the computation speed is greatly improved by using ILU factorization preconditioner for the linear system. In addition, this preconditioning algorithm is simple and parallelizable. Therefore, the algorithm is potential in the applications of weather e-
منابع مشابه
Towards a Cost-effective Ilu Preconditioner with Higher Level Fills
A recently proposed Minimum Discarded Fill (MDF ) ordering (or pivoting) technique is e ective in nding high quality ILU (`) preconditioners, especially for problems arising from unstructured nite element grids. This algorithm can identify anisotropy in complicated physical structures and orders the unknowns in a \preferred" direction. However, the MDF ordering is costly, when ` increases. In t...
متن کاملA Multi-Level Preconditioner with Applications to the Numerical Simulation of Coating Problems
A multi-level preconditioned iterative method based on a multi-level block ILU factoriza-tion preconditioning technique is introduced and is applied to the solution of unstructured sparse linear systems arising from the numerical simulation of coating problems. The coef-cient matrices usually have several rows with zero diagonal values that may cause stability diiculty in standard ILU factoriza...
متن کاملThe flexible incomplete LU preconditioner for large nonsymmetric linear systems
The ILU factorization is one of the most popular preconditioners for the Krylov subspace method, alongside the GMRES. Properties of the preconditioner derived from the ILU factorization are relayed onto the dropping rules. Recently, Zhang et al. [Numer. Linear. Algebra. Appl., Vol. 19, pp. 555–569, 2011] proposed a Flexible incomplete Cholesky (IC) factorization for symmetric linear systems. Th...
متن کاملOn the Scalability of the Albany/FELIX first-order Stokes Approximation ice Sheet Solver for Large-Scale Simulations of the Greenland and Antarctic ice Sheets
We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete ...
متن کاملToward an Effective Sparse Approximate Inverse Preconditioner
Sparse approximate inverse preconditioners have attracted much attention recently, because of their potential usefulness in a parallel environment. In this paper, we explore several performance issues related to effective sparse approximate inverse preconditioners (SAIPs) for the matrices derived from PDEs. Our refinements can significantly improve the quality of existing SAIPs and/or reduce th...
متن کامل